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We propose a simple method for constructing POVMs using any set of matrices which form an 
orthonormal basis for the space of complex matrices. Considering the orthonormal set of irredncible 
spherical tensors, we examine the properties of the construction on the N + 1-dimensional subspace 
of the 2^-dimensional Hilbert space of N qubits comprising the permutationally symmetric states. 
Similar in spirit to Nenmark’s resnlt on realization of a POVM as a projective measnrement, we 
present a method to physically realize the constructed POVMs for symmetric states using the 
Clebsch-Gordan decomposition of the tensor product of irreducible representations of the rotation 
group. We illustrate the proposed construction on a spin-1 system, and show that it is possible to 
generate entangled states from separable ones. 
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Introduction .—Positive Operator-Valued Measures 
(POVMs) are used frequently in quantum information 
processing tasks including quantum filtering [ll 3, 
realization of quantum communication protocols [3l-l6|, 
entanglement verification Q, Remote State Preparation 
(RSP) 0 and conclusive teleportation [l^. Their 
construction and physical realization are of par ticular 
interest in quantum information theory [ll|-[l3j with 
emphasis on optimal measurements through Symmetric 
Informationally-Complete (SIC) POVMs [3, [3| and 
fidelity measures [l6j. Their importance lies in the 
fact that they can be used to achieve certain tasks 
which are outside the scope of projective measurements; 
for example, a set of non-orthogonal states cannot be 
distinguished using projective measurements, but can be 
discriminated unambiguously using POVMs [17|, ll8|. A 
general construction mechanism for POVMs is therefore 
of considerable importance. 

For N qubits residing in the 2^-dimensional Hilbert 
space, the N + 1-dimensional subspace of permutation- 
ally symmetric states corresponds to an experimentally 


developed and rich class 


in quantum information 


processing; the class includes Bell states, Greenberger- 
Horne-Zeilinger (GHZ) states, W states, and V-qubit 
Dicke states, and POVMs on this class are crucial in char¬ 
acterization of nonlocality ^ and multipartite entangle¬ 


ment |22H24j . Given their utility in such tasks, construc¬ 
tion of POVMs whose action always lead to symmetric 
post-measurement states can be particularly useful. 

The purpose then of this Letter is two-fold: to first 
present a simple recipe for construction of POVMs based 
on any orthonormal set of basis matrices on the vector 
space of complex matrices; secondly, use the construction 
to obtain and physically realize POVMs for the permu¬ 
tationally symmetric N -|- 1-dimensional subspace of the 
2^-dimensional Hilbert space in which N qubits reside. 
The space of density matrices and POVMs is a non-linear 
manifold since they are required to be non-negative; how¬ 


ever, it will be seen that the vector space structure of 
the space of complex matrices suffices for constructing 
POVMs. The nature of the construction ensures that the 
POVMs project states to specific subspaces spanned by a 
subset of the orthonormal bases matrices used; the basis 
matrices need not be generators of SU{N), as is usu¬ 
ally assumed in linear representations of the density ma¬ 
trix. Using the representation of the density matrix with 
the orthonormal set of irreducible spherical tensors by 
Fano [ 2 ^, we present a novel method to physically real¬ 
ize the POVMs based on the Glebsch-Gordan decomposi¬ 
tion of tensor products into irreducible representations of 
S't7(2)®^, which include symmetric and anti-symmetric 
states. In the current context, such a procedure is similar 
in spirit to the construction of projective measurements 
from POVMs through Neumark’s theorem. However, the 
POVMs dilated to the 2^-dimensional Hilbert space are 
not projective measurements, and generally not repeat- 
able. Nevertheless, repeatability can be viewed in the 
sense that the dilated POVM always projects onto the 
subspace of symmetric states. We demonstrate the phys¬ 
ical implementation of the POVMs in a 4-dimensional 
2-qubit Hilbert space, and additionally, demonstrate a 
mechanism to generate entangled states from initial sep¬ 
arable states. 

General POVM construction .—We begin with a simple 
observation motivating the construction. Let xi,... ,Xk 
be an orthonormal basis of a fc-dimensional subspace S 
of a vector space V of dimension N, and let X denote the 
N X k matrix whose columns are xi,... ,Xk. Then XX^ 
is Hermitian, non-negative, idempotent, and an orthog¬ 
onal projection operator onto S, spanned by a subset of 
the orthonormal basis xi,...,XNofV. Given any N x N 
matrix A expressed in the basis xi,..., xn, XX^A repre¬ 
sents the projection of A onto S. This presents us with a 
potential POVM construction mechanism for projection 
onto subspaces of the space of density matrices coordina- 
tized by an orthonormal basis. In order to extend this to 
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matrix spaces, we equip the manifold of density matri¬ 
ces with basis matrices for the space ot N x N complex 
matrices; while this ignores the manifold structure of the 
density matrix space, it provides us with a convenient 
mechanism for constructing POVMs. Potential matrix 
bases include, the generalised Gell-Mann basis, Weyl op¬ 
erator basis, spherical tensor basis etc. Such an approach 
for the Bloch vector characterization of the density ma¬ 
trix was adopted in , wherein realization of an entan¬ 
glement witness was considered. 

Let Mn denote the vector space oi N x N matrices 
of dimension over the field of complex numbers C. 
Define a linear map V : Mn ^ Mn, 

V{X) = , 

where Xi is the ith row of X containing N elements, 
and denotes the transpose of a vector v; the map 
transforms an N x N matrix into a column vector of N'^ 
elements by stacking rows into a column, providing an 
equivalent representation of elements oi Mn- The in¬ 
verse map V~^ : Xi IV A4 n is then the map which 
uniquely constructs an N x N complex matrix by tak¬ 
ing the first N elements of V{X) and setting that as 
the first row of X, and so on. Since is a vector 

space over C, we can choose the elementary matrix set 
{Ei^ : z, j = 1,..., N} oi N X N matrices with 1 in the 
(z, j)th position and 0 everywhere else as an orthonormal 
basis; indeed, TrEijE^'^i = and {Eij} forms an 

orthonormal set of basis matrices for XI at . Any element 
oi Mn can be expressed as a linear combination, with 
complex coefficients, of the Eijs. It is easy to verify that 
EijEjj is Hermitian, non-negative, and satisfies 

'£E.jEl = mN. (1) 

i,3 


Let us denote by {T^- : z, j = 1,..., N} another set of 
orthonormal basis matrices for Mn- Now, using the 
map V we can express the orthonormal basis {Eij} as 
{V{Eij) : i,j = 1,..., N}. If U is an N'^ x unitary 
matrix with Ulfi = Ijv 2 , then it is possible to obtain a 
new orthonormal basis {VT^}, where {VT^- = UV{Eij) : 
z,j = l,...,iV}. Explicitly, 


VT.,=C/V(E,,) = (zzg),zzg\ 





(k) 

where the X-dimensional u)/ is the fcth sub-column vec- 
tor of the X^-dimensional complex vector VTij. The in¬ 
verse mapping V“^, along with unitary matrix U, pro¬ 
vides us with way to obtain the N x N orthonormal basis 
matrices {T^} from {Eij} as 


T^J = V-\UV{E,,)) 


( \ 






We have thus shown how any orthonormal basis matri¬ 
ces for M N can be constructed from the elementary basis 
Eij. Note that for any complex matrix A, not necessarily 
square, AA^ is always non-negative and Hermitian. From 
(ED and the fact that UW = Itvs, it is now straightfor¬ 
ward to verify that 


'Eij'I'ij — oIat, 
hi 

where a is a constant not depending on z and j. We 
can formalize the preceding discussion with the following 
Theorem. 

Theorem 1 Let {Tj : j = 1, ..., be an orthonormal 
set of basis matriees for Mn- Then the set {aTjTj} 
satisfies the conditions of a POVM, where a is a constant 
not depending on j. 

It is important to recognize that the subspace onto which 
the POVMs project a state is a vector space and is the 
ambient space in which the submanifold of density ma¬ 
trices reside. 

Spherical tensor representation of density matrix .—A 
density matrix for an fV-qubit system can be represented 
as 


^ (2j + l) ^ ^ ’ 


( 2 ) 


where j = N/2 and are irreducible tensor operators of 
rank k in the V-l-1 dimensional spin space with projection 
q along the axis of quantization in the real 3-dimensional 
space; here Tq =In, the N x N identity operator. Here 
and elsewhere, orthogonality in the matrix space is al¬ 
ways defined with respect to Trace norm or the square of 
the Hilbert-Schmidt norm. The satisfy orthogonality 
and symmetry relations. 


Tr(r,^'rM = (2j + l)<5,,.<5,,,, rf = (-l^r^,. 




where the normalization has been chosen so as to be in 
agreement with Madison convention 27|. The Fano sta¬ 
tistical tensors [1^ or the spherical tensor parameters 
parametrise the density matrix p as expectation values 
of Tg : = tq. In other words, the — 1 spher¬ 

ical tensor operators , in conjunction with the iden¬ 
tity operator, form an orthonormal basis for the vector 
space oi N X N complex matrices over C that acts on the 
V -|- 1 dimensional spin space. In contrast to some or¬ 
thonormal basis matrices, like Gell-Mann matrices, their 
importance lies in the fact they can be constructed as 
symmetrized products of the angular momentum oper¬ 
ators J = {Jx-Jy: Jz) following the well-known Weyl 
construction [2^ as. 


r';{J)=Mkj{J-Vfr^Y^{f), 
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where Afkj are the normalization factors and Y^{f) are 
the spherical harmonics. They possess simple transfor¬ 
mation properties under coordinate rotations of the 3- 
dimensional space: for a rotation R{a,/3,j), where a,/3 
and 7 are Euler angles, the parameters in the rotated 
coordinates, are related to the ones in the initial 

coordinates as 

where DK^{a, [3,^) are the Wigner rotation matrices; 
thus the rank of the tensor is preserved under rotations. 
The matrix elements of the tensor operators, in the an¬ 
gular momentum basis, are given by 

{j'm'\Tg{J)\jm) = [k] C{jkj]mqm'), 

where C{jkj;mqm') are the Clebsch-Gordan coefficients 
and [fc] = The tensor operators are traceless 

but not Hermitian, and cannot in general be identified 
with generators of SU{N). POVMs for permutation- 
ally symmetric states. —The subspace of permutationally 
symmetric states is spanned by the eigen states \jm) of 
angular momentum operators and where j = N/2 
and m = — j,..., -\-j. Employing Theorem [T] in conjunc¬ 
tion with the orthonormal set of spherical tensors from 
@ which form a basis for the space oi N x N complex 
matrices, we can construct POVMs for a symmetric sub¬ 
space of dimension (2j -|- 1) where j = N/2, as 

^ k = 0,...,2j;q = -k,--- ,k, (3) 

where A/” is a constant not depending on k and q. Inter¬ 
estingly, using symmetry properties of Clebsch-Gordan 
coefficients, it can be verified that the set {E^} con¬ 
tains only diagonal matrices. The symmetry relations 
= (—ensure that it is not possible to have 
—1 distinct POVM elements; the POVMs constructed 
are hence neither informationally complete nor symmet¬ 
ric. 

We now present a method which aids the physical real¬ 
ization of the POVM in the laboratory. Recall that N 
qubits reside in a 2^-dimensional Hilbert space of which 
the symmetric space of dimension TV-|-1 is a subspace. In 
this setting, Neumark’s theorem states that any POVM 
on the N I-dimensional symmetric subspace can be 
realized as a a projective measurement in the 2^ dimen¬ 
sional Hilbert space. We adopt an alternative route by 
seeking recourse to Clebsch-Gordan decomposition of ir¬ 
reducible representations of S'f7(2)®^. From the ten¬ 
sor product representation of N qubit state the Clebsch- 
Gordan decomposition leads to a direct sum of irreducible 
representations of which, the one with the largest dimen¬ 
sion corresponds to the invariant subspace of symmetric 
states. In other words, we are able to obtain an orthog¬ 
onal decomposition of the 2 ^-dimensional Hilbert space 


into invariant subspaces including symmetric and anti¬ 
symmetric states, wherein the latter does not transform 
under the action of SU{2)^^. In essence, the decompo¬ 
sition is achieved through the action of a unitary matrix 
with Glebsch-Gordan coefficients as elements. What is 
of interest in this context is that the unitary matrix de¬ 
composes the computational basis in the 2'^-dimensional 
Hilbert space to a set of bases, one amongst which is the 
basis for the N 1-dimensional symmetric subspace. As 
an indirect consequence, we are provided with a mech¬ 
anism to dilate a POVM constructed on the symmet¬ 
ric subspace to the 2'^-dimensional Hilbert space. The 
key difference with Neumark’s theorem is that the di¬ 
lated POVM is not a projective measurement ensuring 
repeatability of measurements; instead, repeatability is 
to be viewed in the sense that measurements based on the 
dilated POVM will always result in a symmetric state. 
Using the constructed POVM, a measurement on an ini¬ 
tial state p® results in a final state p-^ described by the 
density operator. 


Pk,q = Tr{E^p^E^) ■ 

The Spin-1 case .—As an illustration, we now explic¬ 
itly construct the following POVMs for a spin-1 system, 
where the E^s are expressed in the basis |11), |1 0) and 
|1 - 1 ): 



\ 0 0 0 / 


/i 0 0\ /I 0 0 

0 0 0 , Ag = 0 ^ 0 

Vo 0 i/ Vo 0 1 

/ 0 0 0 \ 

Ail = = 0 i 0 , 

Vo 0 iy 
/ 0 0 0 \ 
Aia = 0 0 0 . 

Vo 0 ly 


We can observe the degeneracy regarding the num¬ 
ber of distinct POVM elements; we also note that E^s 
do not transform like spherical tensor operators. We 
now demonstrate how starting with a POVM in a 3- 
dimensional symmetric subspace we can obtain its dila¬ 
tion to the 4-dimensional 2-qubit Hilbert space. Gon- 
sider E/, in the symmetric |lm) basis, m = 1,0,—1. 
The relationship between |lm) basis and the computa¬ 
tional basis is such that |11) = | tt)j 110) = 

11 — 1) = I 4 , 4 ,). Here the spinor in the first and second po¬ 
sitions correspond to the first and second qubits respec¬ 
tively. Let U be the unitary matrix (orthogonal matrix, 
to be accurate, since Clebsch-Gordan coefficients are all 
real-valued) which transforms the computational basis 
to the angular momentum basis |11),|10),|1 — 1), |00). 
Then, the representation of E/ in the 2-qubit state space 
of dimension 4 in the computational basis | tt)j I tDi I it 
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), I U), is given by 


el = U[E\ © 0)C/i, with U = 


/I 0 

0 

7 


1 

©2 

0 0 

0 

1 

75 

1 

©2 

oy 


where © denotes the direct sum. Thus, 


1 

El 


/ 1 0 0 0 \ 
0 i i 0 

0 I I 0 

\ 0 0 0 0 / 


eilV’> = 


/ 1 0 0 0 \ 
0 i i 0 

0 I I 0 

\ 0 0 0 0 / 


/ a \ 
b 

\dj 


and 
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7 (V'kil'0) 7 + 2&c) 


( a \ 

b+c 
2 

b+c 
2 

V 0 y 


where |'0) = a| tt) + ^1 ti) + c| it) + d\ it), with |ap + 
+ |cp + |dp = 1, is the most general pure state in 
the 2-qubit state space. Observe that the normalized 
resultant state is a symmetric state given by 

^(2a2 + 62 + c2 + 26c) ^(2a2 + ^2 + ^2 + 2 bc) 


Therefore, e} projects a vector in the 4-dimensional 
Hilbert space onto the 3-dimensional symmetric space. 
Moreover, in this example, we show that the POVMs can 
be expressed in terms of Pauli spin matrices which are 
identified as Hamiltonians easily implemented in NMR 
quantum computing. If Ii and I 2 represent the identity 
matrices for the qubits 1 and 2 respectively, and cri(l) 
and cri(2), with i = x,y,z, the corresponding Pauli spin 
matrices, then e}’s can be expressed in terms of Pauli 
spin matrices as 


e 


1 

1 


(/l © I2) + © I2) + 2(7 O 0’z(2)) 


where the symbol © denotes the direct product. Similar 
calculations yield: 


1 


Cn = 


^(/l©/ 2 ) + ^( 0 - 71 ) 0 -^ 2 ) 


+ ^(/i © crz(2)) - ^(Ux(l) © o-a:(2)) ; 

^ {h © h) + (ctz( 1) © cr^(2)) 

+ ^K(l) Ocr!/(2)) + ©cra;(2)) ; 


7i = 




g l^(/i © I 2 ) - ^ 771 ) o o 0 - 72 )) 

+ ^7x(l) © (r^{2)) + ^(o-y(l) © <7y{2)) 

^ © h) + ^( 0 - 2 ( 1 ) © 0 - 2 ( 2 )) 

D 


° 612 

5 

12 

1 

24 L 


+ ® + o 7y(l) ® '7y(2)) 


e, = 


--2 


24 L 


{h © I 2 ) + ((^ 2 ( 1 ) © h) 

+ [h © 0 - 2 ( 2 )) + ( 0 - 2 ( 1 ) © 0 - 2 ( 2 )) 
{h © /2) - ((^ 2 ( 1 ) © ^ 2 ) - (/i © ^ 2 ( 2 )) 
+ ( 0 - 2 ( 1 ) © 0 - 2 ( 2 )) + (0-a;(l) © 0-a,(2)) 


We turn our attention to the resulting state following a 
POVM measurement of a spin-1 density matrix . For a 
spin-1 system, the initial density matrix, in the represen¬ 
tation given in (I2|), is given by 

/ 1 + 17 + - 

7 = x -1(7+t 

^ \ 737 

( Pii P12 P13 
P21 P22 P23 
P31 P32 P33 

here k = 1 refers to its vector polarization while k = 2 
refers to its tensor polarization. From (jd]) we obtain. 


L ,2 
7^t0 
2\ 


|(i-i + t-i) 

1 - 727 2 

-1(7+7) 1- 


VstU 

1(^-1 + ^-i) 

17 + ^7 


pfo = 


1 

(pii + P 33 ) 


/ Pll 0 pi3 \ 
0 0 0 
\ P31 0 P33 y 


The only non-zero spherical tensor parameters are 
and t^ 2 - Such a system can be produced in the labora¬ 
tory by the combined electric quadrupole and magnetic 
dipole field, where the direction of the magnetic field is 
along the Z—axis of the Principal Axis frame of the elec¬ 
tric quadrupole field. In similar fashion. 


Pll = 


1 


(pii + P 22 ) 


Pll P 12 0 
P 21 P 22 0 
0 0 0 


Pi-i 


1 

(P22 + P 33 ) 


0 0 0 \ 

0 P 22 P 23 1 ; 

0 P 32 P 33 y 


the non-zero spherical parameters are 77o7±i ^±i- 

In 


(Pll + 16 p 22 + P33) 


( Pll 4pi2 Pl3 
4p21 16p22 4p23 
P31 4p32 P33 
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we can see that all the spherical tensor parameters are 
non-zero. Next, we note that in the case of 

/ 100 \ / 000 \ 

= 0 0 0 , = 0 0 0, (5) 

\ 000 / \ 001 / 


an arbitrary initial state collapses to a separable state 
after performing a measurement with and E^ 2 - 
Entangled states from separable states .— Considering 
a spin -1 example, again, we demonstrate how some of the 
constructed POVM measurements can be used to gener¬ 
ate entangled states from separable states. Consider e} 
acting on a separable state, 


\iP) = 


Then, 


1 

1 

VI/ 


eil^) = 


leading to 


\r) = 


t(l)) + li(l)) 

V2 

/ 1 0 0 0 \ 
^ 0 i i 0 

olio 

V 0 0 0 0 / 

eilV’) _ 1 


t (2)) + I h (2)) 
V2 


V(V'kilV') Vs 


1 

llj 


Cl 

1 

1 

VO/ 


where Itf^) represents the normalized state. In the com¬ 
putational basis, IC-^) can be represented as 

icO = ^[itt) + in) + ut)], 


with corresponding density matrix 


pf — = 


/I 1 1 0\ 
1110 
1110 
V 0 0 0 0 y 


( 6 ) 


In order to check if \tpf) is entangled, we use the PPT 
criterion [s^ , which offers a necessary and sufficient con¬ 
dition for entanglement verification in 2 x 2 and 2x3 
systems. The partially transposed matrix of in ([ 6 ]), 
PpPT^ can be calculated as, 


PppT — 2 


/I 1 1 1\ 
110 0 
10 10 ’ 
\ 1 0 0 0 y 


with eigenvalues Ai = 0.872678, A 2 = —0.333333, A 3 = 
0.333333, and A 4 = 0.127322. Since A 2 is negative we 
conclude that a separable state collapses in to an entan¬ 
gled state after a measurement is performed with e\. 


K.B. acknowledges discussions with Huiling Le about 
Theorem 1. 
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